Theory of Dicke narrowing in coherent population trapping 
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The Doppler effect is one of the dominant broadening mechanisms in thermal vapor spectroscopy. 
For two-photon transitions one would naively expect the Doppler effect to cause a residual broad- 
ening, proportional to the wave- vector difference. In coherent population trapping (CPT), which 
is a two-photon narrow-band phenomenon, such broadening was not observed experimentally. This 
has been commonly attributed to frequent velocity-changing collisions, known to narrow Doppler- 
broadened one-photon absorption lines (Dicke narrowing). Here we show theoretically that such a 
narrowing mechanism indeed exists for CPT resonances. The narrowing factor is the ratio between 
the atom's mean free path and the wavelength associated with the wave-vector difference of the two 
radiation fields. A possible experiment to verify the theory is suggested. 



I. INTRODUCTION 



The spectral line shape of atomic transitions is de- 
termined by many different mechanisms Specifi- 
cally, a Doppler-broadened spectrum can be dramatically 
narrowed due to frequent velocity-changing collisions — 
Dicke narrowing 0, The narrowing factor is pro- 
portional to the ratio between the collisions mean free 
path and the radiation wavelength. Dicke narrowing was 
observed for microwave transitions Q and recently also 
for optical transitions ||. For two-photon transitions, 
such as coherent population trapping (CPT) @, a dra- 
matic narrowing of the expected Doppler- width was also 
observed and was attributed to a Dicke-like narrowing 
B S E3] ■ CPT is a light-matter interaction involv- 
ing thrce-lcvcl atoms and two resonant radiation fields 
(see figure [TJa). The atoms are driven into a coherent 
superposition of the lower levels, |I) and 2) which does 
not absorb the radiation. The efficiency of this process 
strongly depends on the Raman (two-photon) detuning, 
defined as A# = whf — (u>i — u>2)- Therefore, when scan- 
ning the frequency of one of the lasers, an absorption dip 
is observed, which we will refer to as a CPT resonance 
(see figure [I]b). 

The narrow width of CPT resonances is useful for sev- 
eral applications such as frequency standards 0, [Tl| . 
magnetometers [l2|, slow and stored light [H, 13- In 
all the applications a key parameter is the spectral width 
of the CPT resonance. When the CPT is performed in 
room-temperature vapor, the two radiation fields expe- 
rience a Doppler shift which is different for each atomic 
velocity group. In the case of non-degenerate lower lev- 
els (e.g. two hyperfine levels in the ground state of an 
alkali atom) or for non-collinear laser beams, the two ra- 
diation fields experience slightly different Doppler shifts. 
This difference, denoted as the residual Doppler shift, 
results in an effective Raman detuning which is different 
for each velocity group. Therefore, one would naively 
expect a residual Doppler broadening of the CPT reso- 
nance. However, the measured CPT resonance width is 
well below the expected residual Doppler width. For ex- 
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FIG. 1: a. Energy levels scheme of CPT system. Two laser 
fields with frequencies u\ , 0J2 excite two levels to a common 
upper level, b. Typical measurement of a CPT resonance in 
room temperature 87 Rb vapor, with frequency difference of 
lj hf = 6.834GHz. The residual Doppler broadening is about 
lOKHz, while the measured width is only about 300Hz. 



ample, the clock-transition CPT resonance of room tem- 
perature Rubidium vapor is expected to have a residual 
Doppler width of 
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whf « lOKHz, 



(1) 



where i>th is the thermal velocity, whf is the hyperfine 
energy-gap, and c is the speed of light. However, a typical 
measurement of this CPT resonance, depicted in figure 
[T]b, shows a width of only 300Hz. In order to attribute 
this dramatic narrowing to a Dicke-like narrowing effect, 
one commonly assumes that the wavelength which deter- 
mines the narrowing factor is the one associated with the 
frequency difference of the lasers (i.e. the microwave fre- 
quency) p| • In Ref . [IH a numerical model was used to 
calculate the CPT line-shape, introducing discrete veloc- 
ity groups, and indeed demonstrated the expected nar- 
rowing. In the present work we derive an analytic expres- 
sion for the line-shape of CPT resonances, that demon- 
strates both Doppler broadening and collisions-induced 
Dicke narrowing. We show that in a regime relevant to 
most CPT experiments, the narrowing is governed by the 
ratio between the mean free path and the wavelength as- 
sociated with the wave-vector difference of the lasers. In 
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section II we review the two-level Dicke narrowing. In 
section III we derive the theory for Dicke narrowing of 
CPT resonances. Section IV contains a discussion of the 
results and presents a possible experimental setup to ver- 
ify this result. 



II. DICKE NARROWING OF TWO-LEVEL 
ATOMS ABSORPTION LINE 

In this section, as an introduction to the CPT case, we 
develop a theory for the line-shape of a two-level atom, 
including both Doppler broadening and Dicke narrowing. 
Our approach is similar to that presented by Galatry Q . 
We consider a two-level atom with the upper and lower 
states, |3) and |1), and an optical energy gap of Hluq. 
We assume that the kinetic energy of the atom is much 
smaller than hu>o and take the motion of the center of 
mass of the atom to be classical. The atom interacts with 
an external, classical electromagnetic field, with wave- 
vector q and frequency u> = c |q|. We assume that uj is of 
the order of u>o and express the Hamiltonian in the dipole 
approximation and the rotating wave approximation as 



H = hu |3) (3| - hne^-^ |3) (1| + H.c. 



(2) 



where is the Rabi frequency of the field and r = r (t) is 
the time-dependent center of mass position of the atom. 
The equations of motion of the density matrix elements 
p 3i3 and p 3>1 are given by: 

P3,i = -% (uj - iT) p3 A - me^-^ (2p 3 ,3 - 1) (3a) 
p 3 ,3 = -2rp 3 ,3 + 2Imfi*e- i ta- r - wt >p3,i. ( 3b ) 

Here we assumed a single relaxation term (radiation 
bath) , which induces transitions between the atomic lev- 
els (|3) — > |1)) , with the rate T [l6|]. In order to calculate 
the absorption spectrum, we write the energy absorption 
W (i) , in terms of the population response to the applied 
field: 



W (t) = huj a {p 3 , 3 } ficld = 2^o Im (fTe-^ —'Vaa 

(4) 

In the steady-state we take the temporal average, and 
find the absorption spectrum S (u>) to be: 



We consider the unsaturated case in which <C V, 
and use the perturbation expansion of p (t). Taking the 
initial state of the atom to be the ground state, we find 
to zero order in £1 that pfl (t) = pf [ (t) = 0, and to the 
first order that 



with the formal solution 
ft 

p { i\{t)=in dt > e -i{"°-ir){t-t') e i**{t')-i"t\ ( 7a ) 
Jo 

Substituting p^\ (t) in Eq.©, denoting r = t — t' and 
taking the upper limit of the r— integral to infinity, we 
obtain in steady-state, 
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S(w) = Re / dre-^^e^), 
Jo 



(8) 



where A = lj — lu is the detuning from resonance, 

$(r)=q-[r(r)-r(Q)], (9) 

is the phase accumulated during r and the average is over 
atom trajectories, 

e^M = lim - / dtxtaW)-**-^. (10) 

Notice that the accumulated phase includes both the 
Doppler effect and collisions. Assuming that $ (r) is a 
random Gaussian variable and following the cumulant ex- 
pansion procedure [l7j . we substitute e J *( r ) = e~*( T ) 2 / 2 . 
Following Eq.© we write 



Jo 



(11) 



where u a are the velocity cartesian components and a — 
x,y, z, to get 

$(r) 2 = I d< i / d ^ u ° (*i)" Q ' 

Jo Jo 

= 2V g v7 dt (r-i) u a (t) u a ' (0). (12) 
, Jo 

a, a' 

As shown in the appendix, the velocity-velocity corre- 
lation function is given approximately by 



u a {t)u a ' (0) = <W; t 2 h e- 7|t| , 



(13) 



where S a a i is the Kronecker Delta, v t h is the thermal 
velocity and 7 is interpreted as the velocity relaxation 
rate (for the Brownian motion regime) or as the collisions 
rate (for the strong collisions regime). Substituting this 
in Ea.([T2]) gives 



$(r) 2 = 2q 2 vl f dt{r-t)e-^ 
Jo 



^G( 7 r) (14) 



where 

G(x) =x — l + e~ x . (15) 
Then the absorption spectrum, Eq.®, is written as 



3 = -i (c - .T) $ - Xle**'-^ (2p$ - l) , (6) S (u) = J d re -^-r D G( 7 r)/ 7 CQS (At) (lg) 
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FIG. 2: The effect of Dicke narrowing on the absorption spec- 
trum: a numerical solution of Eq. (|16|l for several values of the 
Dicke parameter 2nA/\, with To = 5T. 



where = qv^ = ojv^/c is the Doppler width. Re- 
stricting the discussion to the regime where Tjj ^> T, 
one finds two limits Q: 

(i) The Doppler limit, 7 -C Tp, where the quadratic 
term G (77-) as j 2 t 2 /2 is dominant and the absorption 
spectrum, 



^Doppler (w) 



L ~ 



(17) 



is the well-known Doppler-broadened spectrum. The 
same expression is usually derived by a convolution of 
the homogenous line shape with the thermal velocity dis- 
tribution. 

(ii) The Dicke limit, 7 3> Tjj, where the linear term 
G (77") rs 77" is dominant and the absorption spectrum, 



S 



Dicke 
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A 2 + (r 



r 2 D /7)^ 



(18) 



is a Lorentzian with natural- width T, broadened by 
r|j/7. Denoting the mean free path A = Wth/7, and the 
field's wavelength A = 2ir/q, we find the line- width to 
be r + (27rA/A) Tjj, showing that the Doppler broaden- 
ing is effectively reduced by a factor known as the Dicke 
parameter Q ■ A numerical solution of Eq. (fl6|) , showing 
the transition between the Doppler and the Dicke limits, 
is shown in Fig. ([2]). 



III. DICKE NARROWING OF CPT 
RESONANCES 



A— configuration, with the energy levels E3 > E2 > Ei 
corresponding to an upper state |3) and two lower states 
2} , |1) (see figure [TJa). We will denote the atomic fre- 
quencies uJnn' = (E n — E n i ) I h. The atom interacts with 
external probe and pump electric fields, close to reso- 
nance with the 1 3) — > |1) and |3) — ► |2) transitions, re- 
spectively. The coupling Hamiltonian is 

H c = -hJ2 nne***-^-*"-* |3) (n\ + H.c. (19) 

n=l,2 

where q„ are the wave-vectors of the fields and uj n — 
c|q n |. For brevity, we denote 

Q n (t) = VL n e l ^- r{t ^ w ^ n = 1, 2. (20) 

We consider both dipole relaxation terms, which induce 
transitions between the excited and ground levels |3) — > 
\n) with rates T n (n — 1,2), and within the ground state 
1} <-> |2) with rate r2<->i; and a spin-boson relaxation 
term, which induces adiabatic transitions (decoherence) 
between |1) <-> |2) with rate T ac i- For an atomic vapor 
in the presence of a buffer gas, T n will be the pressure- 
broadened homogenous width. The equations of motion 
for the density matrix are thus 

-2Imft* (t) p 3t i + T 1 p 3t3 + r 2 ^i (p2,2 - Pia) 
-2Imf2 2 (t) p 3i2 + T 2 p 3 ,3 ~ T 2 ^i (p2,2 - Pis) 
2 Im pi (t) p 3 ,i + ^2 (*) P3,2] - (Ti + T 2 ) p 3 ,3 

ift 2 (t) p3,l ~ itt 1 (t) P2.3 - i (W21 - iT 2 l) p2,l 

-iQi (t) (p 3)3 - p lt x) + iVt 2 (t) p2,\ 

i (W3X - iTc) P3,l 

-ifl 2 (t) (p 3 ,3 - P2.2) + itol (t) P\.2 
~ i (^32 - «r C ) P3.2 (21) 

where r c = (r x + r 2 + r 2 ^i + r ad ) /2 , r 21 = iw + 

2T a d- To find the absorption spectrum of the probe field, 
we follow Eqs.([4])-([5]) using the response terms of the 
probe and get, 

T dt p 3<1 (t) 

T fii(i) 

We consider the standard weak probe case, namely 
ill *C il 2 , and avoid saturation by assuming f2 2 smaller 
than any of the relaxation rates r„_^„/ . The perturba- 
tion solution of Eqs. (f2"T|) can easily be done, following 
the same procedure as in the two-level case. To zero or- 
der in ill, with the initial state p(0) = |1) (1| , we find 
p*f\ = 1 and all other matrix elements vanish. To the 
first order in fix, 



Pi.i 
P2.2 

(>3,3 

P\\ 

fa.2 



f 1 

S (u>i) = Im lim / 



(22) 



= iQ 2 (t) p\\[ + iJli (t) - i (w 3 i - *Ti) p\\[ (23a) 
$1 = ifi 3 (t) - i (UJ21 - tT 21 ) p^} (23b) 

and all other matrix elements vanish. We solve formally 
Eq. (|23bjl . assuming the initial state vanishes, 



We now turn to analyze the line-shape of CPT reso- 
nances in the presence of thermal motion and velocity- 
changing collisions. Consider a three-level atom in a 



p^Ut) 



1 j dtiQ,2 (ti) e 
Jo 



(24) 
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and substitute it back into Eq. (|23ap . 



pgi - Mi (t) - i (wsi - iTi) p 3 ]{ (25) 

- ft 2 (t) /" dfeOS (tj) e (-*-«-r»i)(*-*i)4;) , 
Jo 

whose formal solution is an integral equation for (t) , 



(i) 



{1 \ (t) = I tftae (-<WB1_riXt-t!,) [ifii (ta) - 



(26) 



^ 2 (t 2 ) / dtxfiS (*i) e (- i ^- r -)(*=- il ) p ( 1 ) ( tl ) 



An approximate solution for Eq. (|26p is obtained by it- 
erations up to first order in p 3j \. It can be shown, by 
comparing the approximate solution to the exact solu- 
tion for an atom at rest, that the approximation is valid 
in the low-contrast regime, when |fi 2 | < ri 2 r, i.e. m 
the limit of small relative increase in transparency. Note 
that the power-broadening in this regime is small com- 
pared to Ti 2 . 

Performing two iterations of Eq. (|26[) and returning to 
Eq.(|22p. we find the spectrum to be the sum of two terms: 
(i) The one-photon absorption spectrum, 

/>OC 

Si= dre- rir - r ° G ^/T 2 cos(A 1 T) (27) 
Jo 

where Tjj = q^Vth and Ai = u> 3 i — u>i is the one-photon 
detuning of the probe transition. This is essentially the 
two- level absorption spectrum, Eq. (|16[) . and is equivalent 
to the small-signal absorption spectrum in the absence of 
the pump, (ii) The two-photon absorption dip (the CPT 
line-shape) 



f dt r f f 

S 2 = -\n 2 \ 2 Re lim / — / dti / dt 2 / dt 3 
T ^°°Jo 1 Jo Jo Jo 



e (-iAi-ri)(t-ti+t2-t3) e (iA B -r 2 i)(ti-i2) x 
_iq 2 -[r(ti)-r(t 2 )]-iqi-[r(t)-r(i3)] 



(28) 



where An = u>i — lo 2 \ — u 2 is the Raman detuning. For 
brevity, we will take Ai = and write the spectrum in 
terms of A#, i.e. assuming the CPT resonance to occur 
around the center of the one-photon absorption line. The 
calculation can easily be generalized to Ai ^ 0. 

To calculate the shape of the absorption dip, we write 
it as 

S 2 (A fl )=-|fi 2 | 2 Re lim C % fdre^- T ^ 
T ^°°Jo 1 Jo 



dne- TlTl I dr 3 e iK 



(29) 



where 



K = $ 2 (t - t 3 , t) - $1 (t, n + t) 



and $„ (t, t) = q n - [r (t) — r (t — r)] . Both the integra- 
tions over r and t\ contain exponential decay terms, al- 
lowing us to take their upper limits to infinity and aver- 
age the phase-lag term over t, 



poo 

S 2 (A R ) = - \n 2 \ 2 Re / d T e (lA ^ r ^ T 
Jo 



dT ie - riTl I dr 3 e' 



K 



(30) 



In a similar manner to the two-level analysis, we use the 
cumulant expansion to get 



(31) 



Then 



where the first two terms can be found from Eq. (|14| and 
the third term is 



IO = $2 (o, T ) + $2 ( Qj Tl + T ) _ 2$ 2 $! (32) 



$2$! = $ 2 (-t 3 ,t)*i(0,ti+t). (33) 
Using Eqs.(HI]),(IT3l) we find 



*2$i = ^12 / dh / dt 2 u<* (h) u<*' (t 2 ) 



= qi-q 2 u th / dt i 

-Tl — T J — T3 — T 



dt 2 e-^ t2 -^. (34) 



Performing the integral and substituting the results back 
into Eq.([30]). we get 



S 2 (A R ) = - \n 2 \ 2 Re / d re (^-r 21 )T e -(r S 77) 2 G( 7 T) 
Jo 

drie -r 1 r le -(r £ ,/ 7 ) 2 ( 7 r 1 - e -^+e-'-e-'-i) 



f Tl dT3e ni-n 2 v^(e-^-l)(e-^e^+e-^-2)/-y 3 

Jo 

(35) 

where G(x) was defined in Eq. (fT5|) and 

rS s = |qi-q 2 Kh (36) 

is the residual Doppler width. 

Eq. ([33| is a general analytic expression for the CPT 
line-shape. In what follows we limit the discussion to a 
specific realistic regime, which is relevent to most CPT 
experiments. In realizations of CPT in vapor medium, 
the one-photon processes are usually in the far Doppler 
regime (Tfj 3> Ti and Td ^ 7) and in most applications 
the vapor cell contains buffer gas, forcing the two-photon 
processes into the Dicke regime (7 ^ r^ s ^> I\ 2 ). Since 
this regime includes both Doppler and Dicke regimes we 
denote it as the intermediate regime. As described in the 
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previous section, for the Doppler regime we take e _7Tl ~ 
1 — 7n, and thus e -71 " 3 w 1 — 7r 3 and for the Dickc 
regime we take e _7T <C 1 and G (77-) « 7T. With these 
approximations, a simple expression for the CPT line- 
shape is obtained 



S 2 (A*) 



in. 



res 
D 



[r 1 + q r (qi-q 2 )« t 2 h /7] A^- 
where the parameter 



Pi 



17 : 



2tt- 



A 



A: 



CPT 



(37) 



(38) 



is proportional to the ratio between the mean free path, 
A, and the wavelength associated with the wave-vector 
difference, Acpt = 2tt/ |qi — q_2 1 - Equation (|3"T|) is the 
main result of the present work, showing that the CPT 
line-shape is the product of two terms, which are both 
functions of the wave-vector difference. The first term 
determines the line's amplitude and the second is a 
Lorentzian that determines the width. Since 77 multiplies 
the residual Doppler width, it acts as a narrowing factor 
and we denote it as the CPT-Dicke parameter. Finally, in 
a typical setup of a CPT-based frequency standard 
the laser beams are collinear (qi || q 2 ) and the line-shape 
can be written as 



S% (A H ) = 



-m. 



ri2 + ^d s 



(39) 



IV. DISCUSSION AND CONCLUSIONS 

CPT is an inherently narrow-band phenomena usu- 
ally limited by effective broadening mechanisms. Non- 
degenerate CPT resonances in a hot vapor cell with buffer 
gas would naively be broadened by a residual Doppler 
broadening. However, the measured line-width of CPT 
resonances are far below the expected residual Doppler 
width. This effect was attributed to the frequent veloc- 
ity changing collisions with the buffer gas, that are well 
known to narrow atomic absorption transitions in two- 
level atoms. It was also demonstrated using a numerical 
simulation that when frequent collisions occur no Doppler 
broadening is evident in the CPT line [l5j]. In this work 
we developed the theory of Dicke narrowing for CPT res- 
onances in three level atoms. The main result is that the 
residual Doppler width (that would be observable in an 
apparatus with no collisions) is diminished by the ratio of 
the mean free path between collisions and the wavelength 
associated with the wave-vector difference of the two ra- 
diation fields. This theory can be readily extended to 
describe atoms in confined geometries such as thin vapor 
cells and cold atoms traps. 

For hyperfine CPT experiments, performed in vapor 
cells with Alkali atoms and several Torrs of buffer gas, 
the typical CPT-Dicke parameter is 77 ~ 10~ 4 (i.e. very 
strong Dicke narrowing). Hence the residual Doppler 



broadening, which is of the order of a few KHz, is strongly 
reduced and is not measurable (compared to other broad- 
ening mechanisms) . In order to verify the theoretical pre- 
diction it is necessary to increase either the CPT-Dicke 
parameter or the residual Doppler width. An increase of 
the CPT-Dicke parameter can be achieved by decreas- 
ing the effective wavelength or by increasing the mean 
free path between collisions. However, changing the ef- 
fective wavelength is limited by the atomic structure of 
the active atoms and changing the mean free path will 
result in a significant diffusion-broadening. We propose 
to increase the residual Doppler width by introducing a 
small angular deviation, 9, between the pump and probe 
beams. For CPT performed with two degenerate lower 
levels ( I qi I = | q 2 1 ) and small 6, both the residual Doppler 
width and the CPT-Dicke parameter are proportional to 
9. Therefore, the resulting broadening is proportional to 
8 2 and can be increased to a measurable level. 
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APPENDIX 

To obtain the velocity-velocity correlation function of 
Eq. (|I2p we review two simple models: 

(i) In the Brownian motion case the equation of motion 
of the a component of the velocity is 



dt 



u a (t) = - 7 7i a (t) +a a (t) , 



(40) 



where "fu a (t) is the dynamic frictional force on the atom, 
7 is the velocity relaxation rate, and a a (t) is the random 
acceleration [18j . For the correlation function we thus 
obtain 



dt 



u a (t) u a ' (0)= -7it Q (t) u a ' (0) + a a (t) u a (0), (41) 



where the bar indicates ensemble average. Since the ac- 
celeration is not correlated with the initial velocity, the 
last term on the right vanishes, and we have 



u a (t) u a ' (0) = u a (0) u a ' (O)e -71 * 1 



(42) 



In thermal equilibrium with temperature T, the ensemble 
average is taken over the velocity distribution function 



F(u) = {m/2-Kk B Ty v e 



3/2 -mu 2 /2k B T 



(43) 



where m is the mass of the atom, and ks is the Boltz- 
mann factor, and we get 



u a (0) u a ' (0) = 8 a , a v? h , 



(44) 



where 8 a<a i is the Kronecker Delta, and i>th = y/ksT/m 
is the thermal velocity. 

(ii) In the Strong Collisions case the atom is colliding 
with the dilute buffer gas in thermal equilibrium. The 
conditional probability density, f(u,t), to find the atom 
with velocity u at time t, given that at time t — its 
velocity is u (0) , is given by the Boltzmann collision term, 



dt 



dt 



with the initial distribution 

/(u,0) = 5(u-u(0)). 



(45) 



(46) 



The simplest model is the single relaxation rate approx- 
imation, where 



7 (/(M-F(u)), 



(47) 



coll 



F(u) is the equilibrium distribution of Eq . (|4"3")l . and 7 
is the collision relaxation rate. The solution of Eq. (|45|) . 
with Eq. (j47]) . is simply 

/(u, t) = F(u) (l - e-^'l) + e-Tl'l/(u, 0), (48) 



and, with Eq. (|46|) , the velocity of the atom at time t, is 



u(t) = J d 3 uuf(u,t) = u(0)e" 



for consistency 



u a (0)u a ' (0) = / d 3 uu a u a 'F{u) = 



k R T 



(49) 



The velocity- velocity correlation function is 

u a (t) u a ' (0) = u a (0) u a ' (O)e" 71 * 1 . (50) 

and since 

/(u,t-oo) = F(u), (51) 



(52) 



We end up in both these models with 

u a (i) u a ' (0) = WtV" 7 '* 1 . (53) 

where the interpretation of 7 is either the velocity relax- 
ation rate or the collision relaxation rate. 



[1] J. P. Wittke and R. H. Dicke, Phys. Rev. 103, 620 (1956). 

[2] R. H. Dicke, Phys. Rev. 89, 472 (1953). 

[3] L. Galatry, Phys. Rev. 122, 1218 (1961). 

[4] D. Budker, L. Hollberg, D. F. Kimball, J. Kitching, 

S. Pustelny, and V. V. Yashchuk, Physical Review A 

(Atomic, Molecular, and Optical Physics) 71, 012903 

(pages 9) (2005). 
[5] G. Dutier, A. Yarovitski, S. Saltiell, A. Papoyan, D. S. 

ans D. Bloch, and M. Ducloy, Europhysics Letters 63, 35 

(2003). 

[6] E. Arimondo, Progress in Optics, vol. 35 (Elsevier, Am- 
sterdam, 1996). 

[7] N. Cyr, M. Tetu, and M. Breton, IEEE Transactions on 
Instrumentation and Measurement 42, 640 (1993). 

[8] A. Nagel, C. Affolderbach, S. Knappe, and R. Wynands, 
Phys. Rev. A 61, 012504 (1999). 

[9] J. Vanier, M. W. Levine, D. Janssen, and M. Delaney, 
Phys. Rev. A 67, 065801 (pages 4) (2003). 
[10] G. Dutier, P. Todorov, I. Hamdi, I. Maurin, S. Saltiel, 
D. Bloch, and M. Ducloy, Physical Review A (Atomic, 



Molecular, and Optical Physics) 72, 040501 (pages 4) 
(2005). 

[11] S. Knappe, V. Shah, P. D. D. Schwindt, L. Hollberg, 
J. Kitching, L.-A. Liew, and J. Moreland, Applied 
Physics Letters 85, 1460 (2004). 

[12] P. D. D. Schwindt, S. Knappe, V. Shah, L. Hollberg, and 
J. Kitching, Applied Physics Letters 85, 6409 (2004). 

[13] M. D. Lukin, Reviews of Modern Physics 75, 457 (2003). 

[14] M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Re- 
views of Modern Physics 77, 633 (2005). 

[15] M. Erhard and H. Helm, Phys. Rev. A 63, 043813 (2001). 

[16] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Gryn- 
berg, Atom-Photon Interactions (Wiley Interscience, 
New York, 1992). 

[17] M. Abramowitz and I. A. Stegun, eds., Handbook of 
Mathematical Functions with Formulas, Graphs, and 
Mathematical Tables (Dover, New York, 1972). 

[18] S. Chandrasekhar, Astrophys. J. 97, 255 (1943). 



